(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COO PE RATI ON TREATY <PCT) 

llllllllllllllllllllllllllllllll 



(19) World Intellectual Properly 
Organization 

International Bureau 




(43) International Publication Date 
16 December 2004 (16.12.2004) 



PCT 



(10) International Publication Number 

WO 2004/109332 Al 



(51) International Patent Classification 7 : 



G01T 1/166 



(2.1) International Application Number: 

PCT/GB2004/002340 



(22) International Filing Date: 

(25) Filing Language: 

(26) Publication Language: 



3 June 2004 (03.06.2004) 
English 
English 



(30) Priority Data: 

0312776.8 



4 June 2003 (04.06.2003) GB 



(71) Applicant (for all designated States except US): HAM- 
MERSMITH I MANET LTD [GB/GB]; Cyclotron Build- 
ing, Hammersmith Hospital, Du Cane Road, London W12 
ONN (GB). 

(72) Inventors; and 

(75) Inventors/Applicants (for US only): HOGG, Dar- 
ren [GB/AU]; 3/21 Hobart Road, Murrumbeena, VIC 
3163 (AU). THEELEMANS, Kris, Filip, Johan, Jules 
[BE/GB]; Hammersmith Imanet Ltd, Cyclotron Building, 
Hammersmith Hospital, Du Cane Road, London W12 
ONN (GB). 



(74) Agents: HAMMETT, Audrey, Grace, Campbell el al.; 
Amersham plc t Amersham Place, Little Chalfont, Bucking- 
hamshire HP7 9NA (GB). 

(81) Designated States (unless otherwise indicated, for every 
kind of national protection available): AE, AG, AL, AM, 
AT, AU, AZ, BA, BB, BG, BR, BW, BY, BZ, CA S CH, CN, 
CO, CR, CU, CZ, DE, DK, DM, DZ, EC, EE, EG, ES, FI, 
GB, GD, GE, GH, GM, HR, HU, ID, 1L, IN, IS, JP, KE, 
KG, KP, KR, KZ, LC, LK, LR, LS, LT, LU, LV, MA, MD, 
MG, MK, MN, MW, MX, MZ, NA, NI, NO, NZ, OM, PG, 
PH, PL, PT, RO, RU, SC, SD, SE, SG, SK, SL, SY, TJ, TM, 
TN, TR, TT, TZ, UA, UG, US, UZ, VC, VN, YU, ZA, ZM, 
ZW. 

(84) Designated States (unless otherwise indicated, for every 
kind of regional protection available): AR1PO (BW, GH, 
GM, KE, LS, MW, MZ, NA, SD, SL, SZ, TZ, UG, ZM, 
ZW), Eurasian (AM, AZ, BY, KG, KZ, MD, RU, TJ, TM), 
European (AT, BE, BG, CH, CY, CZ, DE, DK, EE, ES, FI, 
FR, GB, GR, HU, IE, IT, LU, MC, NL, PL, PT, RO, SE, SI, 
SK, TR), OAPI (BF, BJ, CF, CG, CL CM, GA, GN, GQ, 
GW, ML, MR, KE, SN, TD, TG). 

[Continued on next page) 



(54) Title: GENERATING DETECTOR EFFICIENCY ESTIMATES FOR A OET SCANNER 




<0 (57) Abstract: Detector efficiency data is generated for a positron emission tomography scanner (2) including a single photon source 
by conducting a blank scan acquisition procedure using the single photon source. The acquired detection count data is processed 
using an efficiency estimation algorithm to calculate data efficiencies of individual detectors in the detector array (8). In one embod- 
iment, the detection count data is output as artificial coincidence count data and the efficiency estimation algorithm operates on the 
artificial coincidence count data. The method can be used in a non-rotating scanner or a rotating scanner. 



WO 2004/109332 Al I Mil 111 111 II II 111 Mill III 1 1! I III II 1 1 III II I II 



Published: 

— with international search report 

— before the expiration of the time limit for amend it ig the 
claims and to be republished in the event of receipt of 
amendments 



For two-letter codes and other abbreviations, refer to the "Guid- 
ance Notes on Codes and Abbreviations" appearing at the begin- 
ning o f each regular issue of the PCI Gazette. 



ill 




WO 2004/109332 



- OEC 2005 



Generating Detector Efficiency Estimates for a PET Scanner 



10 

( 

15 
20 
25 



Field of the Invention 

This invention relates to a method of, and computer software for, 
generating detector efficiency estimates for a positron emission tomography 
(PET) scanner including a detector array for generating detection data and a 
single photon source. 

Background of the Invention 

There are various different types of acquisition procedure in PET, 
including an emission scan, a transmission scan and a blank scan. 

A typical emission scan begins with the injection of a solution including 
a tracer, which is a pharmaceutical compound including a radio-isotope with a 
short half-life, into the subject. The subject may be human or animal. The tracer 
moves to, and is typically taken up, in one or more organs in the subject 
according to biological and biochemical processes which occur within the 
subject. When the radio-isotope decays, it emits a positron, which travels a short 
distance before annihilating with an electron. This annihilation produces two 
high energy photons propagating in opposite directions. The PET scanner 
includes a photon detector array arranged (usually in a ring-shaped pattern) 
around the scanning area. If two photons are detected within a short timing 
window, a so-called "coincidence" is recorded along a line of response (LOR) 
connecting the two detectors. Coincidence counts along each LOR are 
incremented in a data storage part referred to as a sinogram bin each time a 
coincidence is detected between the corresponding detector pair. The output 
coincidence count data from each sinogram bin is typically processed using 
tracer uptake models and image processing techniques to obtain volumetric 
medical images and volumetric tracer uptake rate data for the subject. 

It would be desirable to be able to ascertain detector efficiencies in order 
to be able to detect deterioration in the state of operation of detectors in the 
detector array, since upon deterioration of a detector, the scan results are 
accordingly reduced in accuracy, leading to poor image quality. 
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For quantitative results from PET images, attenuation correction forms 
one of the most crucial data correction stages. In a conventional PET scanner, 
the scanner is provided with one or more positron emitter rod sources, formed of 
a material such as 68 Ge, which emit dual annihilation photons. To derive 

5 attenuation factors, two acquisitions using the rod sources are conventionally 
used - a blank scan, in which the subject being scanned is not present in the 
scanning area (typically, the scanner is empty except for the presence of the 
sources) and a transmission scan in which the subject is present in the scanning 
area. Since the source material is a positron emitter, the two photons arising 

10 from the annihilation of a positron and an electron are acquired in coincidence, in 
the same manner as with an emission scan. Conventionally, the results of the 
blank scan are then divided by the results of the transmission scan to derive an 
attenuation sinogram. The attenuation sinogram is then used to correct the 
emission scan for attenuation, although extra processing steps on blank scan data 

15 and transmission scan data are sometimes used. 

Also, by collecting the coincidence count data during a blank scan and 
processing the count data using an efficiency estimation algorithm, individual 
detector efficiencies can be estimated. For non-rotating PET scanners, estimates 
of the detector efficiency may be calculated by adding together the counts 

20 detected between a particular detector and all the detectors that are in 
coincidence with it. This is known as a fansum and it is based on the following 
measurement model: 

M & = SiSjAij (1) 

where My are the measured coincidence counts between detectors i and j\ 
25 6f and €j are the intrinsic efficiencies of detectors i and j respectively and Ag are 
the ideal coincidence counts measured if the detectors had an ideal performance. 
Thus the measured counts My of a particular coincident pair of detectors are 
proportional to the product of the individual detector efficiencies. An estimate of 
each individual detector efficiency e t is made by summing the measured 
30 coincidence counts My over all detectors j in coincidence with detector i: 
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The above equation is based on the assumption that: 

;A i} —constant (3) 

J 

The assumption is reasonable providing that the efficiencies are not 

significantly different, that£ ^ y i s roughly independent of i and that the number 

J 

5 of coincident detectors is sufficiently large. More accurate methods have been 
proposed which take into account noise and other variables. Further examples of 
efficiency estimation algorithms for calculating detector efficiencies in 
conventional PET scanners are described in "Maximum-Likelihood Estimation 
of Normalisation Factors for PET", D. Hogg, K. Thielemans, T. Spinks and N. 
10 Spyrou, Nuclear Science Symposium and Medical Imaging Conference Record, 
IEEE, 2001. 

In an alternative PET scanner arrangement, used for example in the 
ECAT EXACT3D™ PET scanner, a single photon source is provided for the 
transmission scan. In this scanner the source, a 137 Cs pellet, is automatically 

15 moved by a transport mechanism having a spiral tube that is placed in front of 
the detectors when a transmission scan takes place. By forcing a liquid through 
the tube the source is moved at a constant speed through the scanner in a helical 
motion to provide a full 3D transmission scan. 

As it uses a single photon source however, there is no second photon 

20 generated during an annihilation with which to form a coincident pair. With no 
coincident pair, no real coincidence counts are available from which a 
transmission sinogram could be created. The way this problem is solved is that 
the scanner in this mode measures only the detections furthest from the source, 
ignoring those registered in detectors close to the source, while tracking the 

25 position of the source. The source position is measured at four equidistant points 
about the helix. Since the speed of movement of the source is effectively 
uniform, its position can be calculated accurately to within one detector position 
over the entire length of the helix. From the position of the detector and the 
source, an artificial coincidence count is generated. The set of artificial 

30 coincidence counts detected during a transmission scan are collected and a 



WO 2004/109332 



PCT/GB2004/002340 



4 

transmission sinogram is calculated. However, detector efficiency estimates are 
not obtained from these results. 

Some PET scanners have detector arrangements which rotate during an 
acquisition. Typically, the detectors are arranged in two or more banks which do 
5 not fully surround the subject. Alternatively, the detectors may be arranged in a 
non-ring-shaped pattern. In these cases, the model given by equation (1) above 
is no longer appropriate, as the rotation causes the direct inter-relation between 
detectors and entries in the sinogram bins to break down. Every sinogram bin 
typically contains counts detected by multiple detector-pairs. For some of these 

10 scanners, including the EC AT ART™ PET scanner, the transmission 
measurement is again performed using a single photon point source. However, 
the source is transported only axially relative to the (rotating) detector arrays. 
The EC AT ART scanner contains two such sources, one located on each of the 
opposing detector banks. The ECAT ART scanner is described in further detail 

15 in the article "The ECAT ART Scanner for Positron Emission Tomography: "1. 
Improvements in Performance Characteristics", David W. Townsend et al., 
Clinical Positron Imaging, Vol. 2, No. 1, 1999 and the article "Design and 
Performance of a Single Photon Transmission Measurement for the ECAT ART" 
C.C. Watson, W.F. Jones, T. Brim, K. Baker, K. Vaigneur and J. Young, IEEE 

20 Medical Imaging Conference Record M9-02; 1998. 

A phantom scan involves positioning a body, referred to as the phantom, 
containing a positron-emitting radioactive source material in the scanner. A 
phantom scan is either an emission or a transmission scan with a 'phantom* in the 
scanner instead of a patient. Typically, the 'phantom* is an object made out of 

25 plexiglass or suchlike and filled with water mixed with a radioactive substance. 
The phantom has known shape and attenuation characteristics. A phantom 
emission scan has in the past been used to calculate detector efficiency estimates. 
A problem with having to conduct a phantom scan is that the phantom must be 
handled by an operator to place the phantom inside the scanner and to 

30 subsequently remove the phantom. It is both inconvenient and time-consuming, 
since in any event an operator needs to be available to handle the phantom at the 
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start and end of the procedure. It is also potentially hazardous, in terms of lifting 
the phantom and exposure to radioactivity. 

It would be desirable to be able to estimate detector efficiency data in 
PET scanners of a type which use a single photon source for transmission scans, 
5 without the drawbacks of using a phantom scan. 

In accordance with one aspect of the present invention there is provided a 
method of generating detector efficiency data for a positron emission 
tomography scanner including: 

a detector array for generating detection data; and 
1 0 a single photon source, 

wherein the method comprises: 

conducting an acquisition procedure using the single photon source to 
produce detection data; and 

processing said detection data using an efficiency estimation algorithm to 
15 calculate data representative of the efficiencies of individual detectors in said 
array. 

By enabling the estimation of detector efficiencies from an acquisition 
using the single photon source, detector efficiencies can be generated without 
significant inconvenience to an operator. The detector efficiencies may be 
20 derived from a blank scan acquisition conducted at the operator's convenience. 
Furthermore, the regular need for the use of a phantom scan procedure can be 
avoided. 

In a preferred embodiment of the invention, detector efficiency estimates 
are made using the artificial coincidence counts generated during a blank scan 

25 acquisition made using the single photon source. Typically, artificial 
coincidence counts are the only suitable detection data made available as an 
output from a scanner of a type such as an ECAT EXACT3D PET scanner. 
However, detector efficiency data cannot be accurately estimated from artificial 
coincidence counts using known techniques, because the known measurement 

30 models do not apply. The present invention provides a new measurement model 
and exemplary efficiency estimation algorithms, which can be applied to 
artificial coincidence counts produced using blank scans. 
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Further features and advantages of the present invention will become 
apparent from the following description of preferred embodiments of the present 
invention, made by way of example only with reference to the accompanying 
drawings. 

5 

Brief Description of Drawings 

Figure 1 is a schematic diagram of a PET scanning facility; 

Figure 2 is a schematic diagram of data processing components in a non- 
rotating PET scanner and an associated data processing terminal; 
10 Figure 3 is a schematic cross section of the arrangement within a non- 

rotating PET scanner during an emission scan acquisition; 

Figure 4 is a schematic cross section of the arrangement within a non- 
rotating PET scanner during a blank scan acquisition using a single photon 
source; 

15 Figure 5 is an illustration of a user interface provided in an embodiment 

of the invention for presenting detector efficiency estimates to an operator; and 

Figure 6 is a schematic cross section of the arrangement within a rotating 
PET scanner during a blank scan acquisition using a single photon source. 

20 Detailed Description of Preferred Embodiments 

Referring now to Figure 1, a PET scanning facility, arranged in 

accordance with an embodiment of the invention, is shown, which includes a 

PET scanner 2, a tracer generator module 4 and an operator computer terminal 6. 

The scanner 2 includes a detector array 8 arranged about a scanning area, in 
25 which a subject 10 is located during a transmission scan and during an emission 

scan. 

As shown in Figure 2, the PET scanner 2 includes a control unit 30, 
detection data processing circuitry 32, one or more single photon source position 
detectors 34, a count store 36 and an Input/Output (I/O) device 38. The 
30 computer terminal 6 includes a central processing unit (CPU) 42, memory 44, 
hard disc drive 46 and I/O device 40, which facilitates interconnection of the 
computer 6 with the PET scanner 2. Operating system programs 48 are stored on 
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the hard disc drive 46, and control, in a known maimer, low level operation of 
the computer terminal 6. Program files and data 50 are also stored on the hard 
disc drive 46, and control, in a known manner, outputs to an operator via 
associated devices. The associated devices include a display, a pointing device 
5 and keyboard (not shown), which receive input from, and output information to, 
the operator via further I/O devices (not shown). Included in the program files 
50 stored on the hard drive 46 are a detector efficiency calculating application 52 
and a detector efficiency user interface application 54. A database 56 is used to 
store the detection data transferred from the PET scanner 2. 

10 The detection data processing circuitry 32 of the PET scanner operates in 

two different modes. In a positron source mode, the detection data processing 
circuitry 32 processes all events detected in the detector array 8, and by using a 
coincidence timing window, detects coincidences between the events which are 
recorded as real coincidence counts My. These real coincidence counts are then 

15 output to sinogram bins in the count data store 36 for subsequent processing. 

In a single photon source mode, the detection data processing circuitry 32 
generates artificial coincidence counts when a single photon source is exposed in 
the scanner, as will be described below in further detail. The one or more single 
photon source position detectors 34 provide a position detection output to control 

20 unit 30 which generates a dynamic position estimate during the acquisition, 
which is sent to the detection data processing circuitry 32. The detection data 
processing circuitry 32 filters the detector outputs so as to discard all events 
detected on the side of the detector array nearest the current source position. 
Next, the detection data processing circuitry 32 generates an artificial 

25 coincidence count M'g 9 and transmits the artificial coincidence count data to 
sinogram bins in the count data store 36 for subsequent processing. 

One implementation of the invention, which is described in further detail 
below, relates to a non-rotating PET scanner. In the embodiment of this 
implementation, a PET scanner of the ECAT EXACT3D™ type is used. 

30 However, it should be noted that this implementation of the invention applies 
also to other non-rotating scanners which include a single photon source. In this 
embodiment, the detectors in the detector array 8 are arranged in square detector 
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blocks, each containing 64 detector elements. The detectors are arranged in six 
rings of detector blocks. There are 48 detector blocks in each ring. The total 
number of lines of response (LORs) that can be acquired is thus, in this 
embodiment, about 190 million. 
5 Figure 3 illustrates features of operation of the non-rotating PET scanner 

in a positron source mode during an emission scan. During an emission scan, the 
subject 10 is placed in the scanning area and contains the tracer generated by the 
tracer generator module 4. Figure 3 shows a positron emission event being 
registered in the detector array 8. The positron 60 annihilates and generates a 

10 first photon 62 travelling in one direction and a second photon 64 travelling in an 
opposite direction. The first photon is detected by a detector element 70 on one 
side of the detector array 8, and the second photon 64 is detected in a different 
detector element 72 on the other side of the detector array 8. As described 
above, the detection data processing circuitry 32 registers the two as a 

15 coincidence along the LOR defined between the two different detectors 70, 72. 

Figure 4 illustrates features of operation of the non-rotating PET scanner 
in a single photon source mode during a blank scan. In this mode, a single 
photon source 80 is transported around the periphery of the scanning area using a 
single photon source transport mechanism 82. The single photon source 80, 

20 referred to as a "point" source due to its small diameter, consisting of a pellet of 
Cs is driven in a fluid-filled steel tube wound into a helix and positioned just 
inside the detector array 8. The tube diameter is 4mm or smaller and the 137 Cs 
point source, with a radioactivity typically between 1 to 10 mCi, has a diameter 
of 2mm or less. The source is transported in a liquid carrier and is buoyancy 

25 matched to the specific gravity of the liquid. The arrangement is designed to 
increase consistency of the speed of travel of the source. 

There are 48 turns of the helix within the detector field of view, 
coiresponding to one turn per detector ring. The source carrier is driven by a 
pump (not shown) at a speed of about 1ms" 1 . The source 80 is shown in Figure 4 

30 emitting a single photon 90 which is detected at a detector 94 on the opposite 
side of the detector array 8 to the source 80. Figure 4 shows a virtual second 
photon 92, representing the artificial coincidence registered in the detector data 



WO 2004/109332 



PCT/GB2004/U02340 



9 

processing circuitry 32 as described above. An event along an LOR is registered 
between the two respective detectors involved, the detector 94 registering the 
event and the detector 96 immediately adjacent the current source position 80. 
As described above, the detection data processing circuitry 32 generates an 
5 artificial coincidence count which is transmitted to the count data store 36. 

After the end of the acquisition, the operator can instruct computer 
terminal 6 to download the detection data held in count data store 36 to the 
database 56 on the computer terminal hard drive 46, where the data can be 
processed. 

10 hi accordance with this embodiment of the present invention a 

measurement model for an acquisition using the single photon source is 
provided, which takes into account the method whereby the artificial coincidence 
counts are produced. The measurement model has the form: 

M'^SiA'ij+EjA'ji (4) 

15 where M'y are the measured artificial coincidence counts between detectors i and 
j 9 6/ and €j are the intrinsic efficiencies of detectors i and j respectively and A 
are the ideal artificial coincidence counts measured if the detector efficiency of 
detector i had an ideal performance, while the detection efficiency of detector j 
was 0. This measurement model takes into account that an event M ' £J recorded 

20 between detectors i and j will have originated either from detector i or j. These 
will have each have a respective inherent efficiency or 6j. For a ring-shaped 
scanner where the single photon source orbits a spiral with a centre of rotation 
coinciding with the axis of the PET scanner, it can be assumed that A )j= A y f . 

Thus, the measurement model used is additive, in that the measured 

25 artificial coincidence counts M'y of a particular artificially coincident pair of 
detectors are equal to a weighted sum of their individual efficiencies. This 
contrasts with the known model for a true coincidence system, based on equation 
(2) above, which is multiplicative, in that the measured counts My of a particular 
coincident pair of detectors are proportional to the product of the individual 

30 detector efficiencies, which can lead to very significantly different results. 
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Next, an algorithm is provided for estimating detector efficiencies. 
Assuming that the transmission scan is reasonably constant over the region of the 
coincident pairs taken into consideration, then the A'g terms can be ignored. 
Equation (4) can then be rearranged and solved as an efficiency estimation 
5 algorithm using the following iterative formula: 
( \ 

S^=± 2>V-Z4 k=l,.^,Nuer (5) 

V J J J 

where J is the number of detectors in the accepted fan and k is the iteration 
number. For the first iteration the efficiencies can be initialized to any positive 
constant. For a fast convergence to a solution, new estimates of sf +l are passed 

10 into the right hand side of equation (5) as soon as they have been calculated. If 
the A *ij terms cannot be ignored, a general iterative linear equation solver may be 
used. It is possible to use more sophisticated methods, for example based on a 
Maximum Likelihood prescription which takes noise in the measurements into 
account. In particular, because of the additive nature of the measurement model 

15 (4), the Maximum Likelihood methods applicable to image reconstruction in 
PET, such as the MLEM algorithm as described in the article "Maximum 
likelihood reconstruction for emission tomography", Y.Vardi, L.A Shepp, IEEE 
Trans. Med. Imag., Vol. 1, pp.113-121,1982, may be used. To briefly describe 
the MLEM algorithm, we introduce some notation as follows. We use the index 

20 z to run over all sinograms bins. We define a matrix P to concisely write the 
linear measurement model as: 

Mz^P^j (6) 
J 

or in matrix notation: 

M -Ps (7) 
25 The MLEM algorithm can then be written as: 

s k+\ +Pu)®P t \p[+Pe k ) (8) 

where u is a sinogram with all elements equal to 1, + denotes element-wise 
division and 0 denotes element-wise multiplication. The initial estimate for s is 
arbitrary, as long as all elements are strictly positive. However, we recommend 
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setting all s j to 1 . For more details, we make reference to the above mentioned 
article, the contents of which are incorporated herein by reference. 

Whichever efficiency calculating algorithm is selected, the detector 
efficiency calculating component 52 implements the efficiency calculating 
5 algorithm as provided hereby. 

A PET scanner can perform a transmission scan in either 2D or 3D mode. 
For patient scans on the ECAT EXACT 3D, both 3D blank and 3D transmission 
acquisitions can be conducted. Data from several LORs are then combined 
together in a process called spanning to save data storage space. This means that 

10 the above measurement model (4) is no longer appropriate. To avoid this 
problem, a blank scan acquisition is preferably conducted in a 2D mode (in 
which no spanning occurs), and from this individual detector efficiencies can be 
calculated. However, where it is possible to acquire the blank scan acquisition 
data in 3D mode without spanning the data, the 3D blank scan data may be used 

15 to calculate detector efficiencies in accordance with the invention. Even with 3D 
blank scan data with spanning, although more complicated, the measurement 
model is still a linear combination of the detector efficiencies, and hence this 
system of linear equations can be solved by an iterative procedure similar to 
equation (5) above. Similarly, a Maximum Likelihood approach such as the 

20 MLEM algorithm can be applied. 

To experimentally test the efficiency estimate algorithm provided using 
equation (5) above, a direct plane 1 hour blank scan acquisition was performed 
on an ECAT EXACT3D scanner. The efficiencies were calculated using 
equation (5) and a corrective efficiency sinogram generated for each ring. The 

25 acquired blank sinogram, generated using the artificial coincidence count data, 
showed clearly the differing detector efficiencies by diagonal lines of varying 
intensity. The acquired blank sinogram was corrected for these effects using the 
calculated efficiency sinogram, and the corrected blank scan sinogram was seen 
to be free of diagonal lines of varying intensity. Thus, the calculated efficiency 

30 estimates were shown to be accurate. 
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Figure 5 illustrates a display part of a user interface generated on the 
computer terminal 6 by the detector efficiency user interface application 54 using 
the detector efficiency data calculated by data efficiency calculating application 
52 after a blank scan procedure as described above. The display includes a 
detector representation 100, showing each of the detectors in the detector array 8 
on a flat grid 100. The flat grid 100 is divided into detector block squares 102, 
representing each detector within the block as a separate pixel for the different 
brightness intensity. The brightness intensity used in each pixel represents the 
calculated efficiency estimate for the respective detector. A relatively dark pixel 
104 is used to illustrate a detector having a relatively poor efficiency. A 
relatively bright pixel 106 is used to illustrate a detector having a good, folly 
operative efficiency. An intermediate brightness intensity pixel 10S is used to 
represent a detector having an intermediate efficiency. Thus, by viewing the grid 
100, the operator is able to visually identify detectors having poor and 
intermediate efficiencies, thereby to detect problems in the detector array 8. For 
example, occasionally an entire detector block will be somehow impaired or 
even inoperative. In this case, the entire detector block square 102 will appear 
relatively dark, and the operator can determine that the block requires 
replacement in order to keep the PET scanner image quality high. 

Furthermore, the detector efficiency user interface application 54 is 
provided with one or more automated block impairment check routines which 
are used to highlight blocks which are determined to be below a given threshold 
of acceptable operation. In a first check routine, an average efficiency is taken 
across each block. If the average detector efficiency within any block is below a 
preset threshold, the relevant block square 102 is highlighted, for example by the 
addition of a coloured border or a flashing display of the block square 102. 

In a second block check routine, the coefficient of variation of 
efficiencies within a block is calculated. The detector efficiency user interface 
application 54 then highlights, for example by the addition of a coloured border 
or a flashing display, any detector square 102 representing a block having a 
coefficient of variation which is higher than a preset threshold. 
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Other implementations of the invention, which are described in further 
detail below, relate to a rotating PET scanner. In an embodiment of this 
implementation, a scanner of the ECAT ART™ type is used. However, it should 
be noted that this implementation of the invention applies also to other rotating 
5 scanners which include a single photon source, hi this particular embodiment, 
the detectors in the detector array are arranged in square detector blocks, each 
containing 64 detector elements. The detectors are arranged in two banks of 
detector blocks. The banks do not entirely surround the scanning area, thus 
allowing fewer detector elements to be used in the scanner. 

10 Figure 6 illustrates features of operation of the rotating PET scanner in a 

single photon source mode during a blank scan. The scanner has first and second 
detector banks 120, 122 arranged on different sides of the patient. The two 
detector banks are arranged on a rotating mechanism which allows the banks to 
rotate about the scanning area. For blank scans, a single photon source 124, 130 

15 is held in a single photon source transport mechanism 126, 132 on one side of 
each of the respective detector banks 120, 122. The single photon sources 124, 
130 are thus mounted to, and rotate with, the detector array. Each single photon 
source transport mechanism is arranged to transport the source axially of the 
rotating detector array during a blank scan. A single photon source shield 128, 

20 134 is provided for each single photon source. During a blank scan, one or each 
of the single photon source shields 128, 134 may be retracted to expose the 
single photon source in the interior of the scanner. Figure 6 shows an 
arrangement in which the first single photon source 124, on the first detector 
bank 120, is exposed by the retraction of the first single photon source shield 

25 128, and in which the second single photon source 130, on the second detector 
bank 122, is exposed by the retraction of the second single photon source shield 
130. In this arrangement, a photon 138 emitted by the first single photon source 
124 is detected at a detector element 140 on the second detector bank 122. 

For a rotating scanner, a count in a given sinogram bin corresponding to a 

30 single LOR, may correspond to different detector pairs, because when the 
detector array is rotating, the LOR between two detectors at time 0 would later 
be covered by another detector pair. 
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In one embodiment of the invention, detector efficiency estimates are 
obtained by conducting a blank scan acquisition when the rotatable detector 
array is held stationery. In the case that no LORs are combined (i.e. no 
spanning), there is a one-to-many correspondence between the (non-zero) bins in 
5 the sinogram and the detectors. In this case, the measurement model is: 

M\=e h A\ (9) 

where there is a one-to-many correspondence between j z and z, and 
efficiencies can be solved by using: 

10 where the set Z y contains all sinogram bins that correspond to the 

detector/, te. 

Zj={z\j z =j} (11) 

Equation (10) is the Maximum Likelihood solution for the efficiencies 
when the measurement model (9) applies. The sum in (10) over all sinogram 
15 bins that correspond to the same detector reduces the noise in the estimated 
efficiencies. 

In further embodiments of implementation of the invention for rotating 
scanners, detector efficiency estimates are generated when the detector array is 

20 rotating. The detector efficiencies may depend on temperature, and hence 
rotation speed, thus detector efficiency estimates may be generated which more 
accurately reflect the operative efficiencies. 

In one embodiment in which the detector array is rotated during a blank 
scan acquisition, one single photon source is exposed to the (rotating) detector 

25 array, and the other single photon source is shielded, as shown in Figure 6. In 
this case, a measurement model that can be used has the same fomi as equation 
(9) above, albeit with different A\ values, where there again is a one-to-many 
correspondence between j z and z. Hence a prescription similar to (10) can also 
be used. The measurement is then to be repeated with the other source viewable 

30 such that the other detector bank is irradiated. 
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In a further embodiment in wliich the detector array is rotated during a 
blank scan acquisition, both single photon sources are exposed in the (rotating) 
detector array. In this embodiment, at least some of the sinogram bin values will 
correspond to either two or four overlapping but different source-detector 

5 arrangements. A possible solution to this problem is to use gating hardware. 
Gating hardware provides multiple sinograms (one for each gate), where counts 
are assigned to a gate according to the time they are detected. The overlap effect 
occurs at the periphery of the detector banks. There is an angle, 5, which defines 
the smallest rotation angle before overlap occurs. Taking the largest angle 5' <6 

10 such that 360° is an integer multiple of d\ i.e. 360° = 5' * TV, where N is an 
integer, N gates are set up, each of length T/N where Tis the period of rotation of 
the scanner. After a full rotation, the scanner hardware will add counts to the 
same set of gates again, such that N different sinograms, each with uniquely 
identifiable correspondence between the LOR and detector, are produced. A 

15 measurement model, and analysis function, such as that described above using a 
single exposed source may then be employed to generated detector efficiencies. 

Finally, also for the case of rotating scanners, spanning is often 
introduced to reduce storage demands. As for non-rotating scanners, this will 
mean that the measurement model (9) is to be modified such that the measured 

20 artificial counts are a linear combination of efficiencies. As before, one can solve 
for the efficiencies with a Maximum Likelihood algorithm such as the MLEM 
algorithm. 

Thus, it will be appreciated that the present invention may be 
25 implemented to generate detector efficiencies from a blank scan that could then 
be used as a daily check on the performance of a PET scanner. The invention 
may also be used to provide detector efficiency factors for a scanner 
normalisation procedure. Since detector efficiencies change over time, it is 
desirable to be able to ascertain detector efficiency factors on a regular basis and 
30 without the need for complex procedures. The above-mentioned efficiency 
estimation algorithms, or alternatives thereto, may be used in such a procedure. 
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The above embodiments are to be understood as illustrative examples of 
the invention. It is to be understood that any feature described in relation to any 
one embodiment may be used alone, or in combination with other features 
described, and may also be used in combination with one or more features of any 
5 other of the embodiments, or any combination of any other of the embodiments. 
Furthermore, equivalents and modifications not described above may also be 
employed without departing from the scope of the invention, which is defined in 
the accompanying claims. 
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Claims 

1. A method of generating detector efficiency data for a positron 
emission tomography scanner including: 

5 a detector array for generating detection data; and 

a single photon source, 
wherein the method comprises: 

conducting an acquisition procedure using the single photon source to 
produce detection data; and 
10 processing said detection data using an efficiency estimation algorithm to 

calculate data representative of the efficiencies of individual detectors in said 
array. 

2. A method according to claim 1, where said acquisition procedure 
1 5 includes a blank scan acquisition. 

3. A method according to claim 1 or 2, wherein the scanner includes 
a coincidence detection system for producing coincidence count data (My) in the 
detection data during an acquisition when a positron source is inside the scanner, 

20 and wherein the scanner is arranged to produce artificial coincidence count data 
(M'iJ) during an acquisition using the single photon source, 

and wherein the step of processing said detection data comprises 
processing said artificial coincidence count data. 

25 4 - A method according to claim 3, wherein the efficiency estimation 

algorithm is based upon a measurement model which is additive, in that the 
measured counts of a particular artificially coincident pair of detectors is related 
to a weighted sum of their individual efficiencies. 

30 5 - A method according to any preceding claim, wherein the scanner 

is a non-rotating scanner. 
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6. A method according to any of claims 1 to 4, wherein the scanner 
is a rotating scanner. 

7. A method according to claim 6, wherein the scanner comprises 
5 two single photon sources and the method comprises selectively operating one of 

the two single photon sources during the acquisition procedure. 

8. A method according to any preceding claim, further comprising 
generating an output, responsive to said data representative of efficiencies, on an 

10 output device for an operator. 

9. A method according to any preceding claim, comprising 
processing said data representative of efficiencies to identify detector elements, 
or groups of detector elements having relatively low efficiencies. 

15 

10. A method according to claim 9, comprising processing said data 
representative of efficiencies using a function determining a parameter relating to 
an average over a plurality of detector elements. 

20 11. A method according to claim 9 or 10, comprising processing said 

data representative of efficiencies using a function determining a parameter 
relating to an amount of variation therein. 

12. Computer software for generating detector efficiency data for a 
25 positron emission tomography scanner including: 

a detector array for generating detection data; and 
a single photon source, 

wherein the software is adapted to operate on detection data generated by 
conducting an acquisition procedure using the single photon source, and 
30 wherein the software is adapted to process said detection data using an 

efficiency estimation algorithm to calculate data representative of the efficiencies 
of individual detectors in said array. 
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13. Computer software according to claim 12, where said acquisition 
procedure includes a blank scan acquisition. 

5 14. Computer software according to claim 12 or 13, wherein the 

scanner includes a coincidence detection system for producing coincidence count 
data (My) in the detection data during an acquisition when a positron emitting 
source is inside the seamier, and wherein the scanner is arranged to produce 
artificial coincidence count data (M'ij) during an acquisition using the single 

10 photon source, and wherein the software is adapted to operate on said artificial 
coincidence count data. 

15. Computer software according to claim 14, wherein the efficiency 
estimation algorithm is based upon a measurement model which is additive, in 

15 that an efficiency of a particular artificially coincident pair of detectors is related 
to a sum of their individual efficiencies. 

16. Computer software according to any of claims 12 to 15, wherein 
the scanner is a non-rotating scanner. 

20 

17. Computer software according to any of claims 12 to 15, wherein 
the scanner is a rotating scanner. 

18. Computer software according to claim 17, wherein the scanner 
25 comprises two single photon sources and the method comprises selectively 

operating one of the two single photon sources during the acquisition procedure. 

19. Computer software according to any of claims 12 to 18, wherein 
the software is adapted to generate an output, responsive to said data 

30 representative of efficiencies, on an output device for an operator. 
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20. Computer software according to any of claims 12 to 19, wherein 
the software is adapted to process said data representative of efficiencies to 
identify detector elements, or groups of detector elements having relatively low 
efficiencies. 

5 

21. Computer software according to claim 20, wherein the software is 
adapted to process said data representative of efficiencies using a function 
determining a parameter relating to an average over a plurality of detector 
elements. 

10 

22. Computer software according to claim 20 or 21, wherein the 
software is adapted to process said data representative of efficiencies using a 
function determining a parameter relating to an amount of variation therein. 

15 23. A data carrier comprising computer software according to any of 

claims 12 to 22. 
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